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Abstract 

Inspired by dengue and yellow fever epidemics, we investigated the contact process (CP) in a 
multiscale network constituted by one-dimensional chains connected through a Barabasi-Albert 
scale-free network. In addition to the CP dynamics inside the chains, the exchange of individuals 
between connected chains (travels) occurs at a constant rate. A finite epidemic threshold and 
an epidemic mean lifetime diverging exponentially in the subcritical phase, concomitantly with 
a power law divergence of the outbreak's duration, were found. A generalized scaling function 
involving both regular and SF components was proposed for the quasistationary analysis and the 
associated critical exponents determined, demonstrating that the CP on this hybrid network and 
nonvanishing travel rates establishes a new universality class. 
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I. INTRODUCTION 



Complex networks are involved in a wide range of biological, social and technological 
systemshj, 2| such as, for instance, the Internet ll, 3|, the network of human sexual contacts 

fin n n 

[4j,|5j, and the transportation infrastructure |6|,|7j]. These nets exhibit "scale-free" (SF) degree 
distributions, characterized by a probability P{k) ~ k~ J that an element in the network is 
connected to k other elements. Usually, the degree exponent assumes a value in the range 
2 < 7 < 3 In special, since many infectious diseases exploit the social contact network 
- o., t0 ^ ^ B - — r ^e S ^ 8ignlficant damages 

throughout the Internet [9j, epidemic spreading through a variety of networks, including 
random graphs (RG), small worlds (SW) and SF networks, have been extensively studied 
9rll5|. 

Pastor-Satorras and coworkers investigated the epidemic spreading on complex 

networks using the susceptible-infected-susceptible (SIS) model. For all networks considered 
(RG, SW, and SF), excellent agreements between mean-field and numerical results were ob- 



tained at the transition to the absorbing state 
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111 ]. However, Castellano and Pastor-Satorras 



14i | observed a non-mean-field critical behavior dependent on 7 for the classical Harris con- 
tact process (CP) 16j| on SF networks. So, these models, sharing the directed percolation 
(DP) universality class on regular lattices j3, exhibit distinct critical behaviors on SF 
networks. 

The deep impact of complex networks on epidemiology, radically altering traditional 
control and vaccination strategies, bring to the focus some epidemic processes associated 
to tropical diseases that can not be properly described by either a purely regular, SW 
or SF networks. Indeed, the spreading of dengue and yellow fever, the main re-emergent 
diseases [18] in the world, are determined by the population dynamics of mosquitoes with 
a limited flight range (~ 200 m for Aedes aegypti), domestic habits and local dispersion 
independent on the human contact network. Furthermore, in their epidemiological chains, 
a relevant factor for viral circulation is the travels of symptomless infected people. Such 
travels occur on regional and national interurban commuting traffic networks which behave 
as SF weighted nets with modularity (7]. Hence, the spreading of dengue and yellow fever 
combines a local contact process among mosquitoes and individuals around their houses with 
long-range dispersal through symptomless infected people traveling between cities. Aiming 
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to include this two-component spreading process, Silva et al. studied the SIS model in SF 



networks of square lattices 



191 ] . Although they focused on dynamics instead of the transition 



to the absorbing state, a nonzero epidemic threshold characterized by a singular approach 
to the epidemic-free state was found. 




FIG. 1: Schematic representation of the network. The connections between chains (the cities) are 
represented by lines. Black and white squares correspond to occupied and empty sites, respectively. 



In this paper, we report on numerical simulations of the CP model in a SF network in 
which the nodes are linear chains. The transition to the absorbing state is focused. The 
model and its computer implementation are presented in Sec. [Ill Simulations are reported 
and discussed in Sec. IHIi Finaly, some conclusions are drawn in Sec. IIVI 



II. MODEL 



In the model, individuals lie on iV linear chains of size L with periodic boundary conditions 
(the cities), in which empty and occupied sites represent healthy and infected individuals, 
respectively. The cities are connected according to the Barabasi- Albert (BA) model lj], i. 
e., in a SF network with a degree distribution given by P(k) ~ k~ 3 (FIG. [1]). The model 
dynamics incorporates a local spreading in which healthy sites with n occupied nearest 
neighbors (NNs) are infected at a rate n\/2 while any infected site is spontaneously cured 
at rate 1. Additionally, two sites in distinct connected nodes are interchanged at a rate a 
representing a travel. These travels occur preferentially to nodes with larger connectivity. 
The simulations were implemented as follows. At each time step, an infected site is selected 
at random and the time incremented by At = l/N occ , where N occ is the total number of 
infected sites in all nodes. The selected site will perform one of three actions: (i) become 
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healthy with probability p = 1/(1 + A + a); (ii) randomly select and infect one of its healthy 
NNs with probability q = + \ + a); (iii) travel with probability r = 1 — p — q. In a 
travel, the target node j, with degree kj, is chosen among all those connected to the departure 
node i with probability IL-^ = kj/Y2j kj- Notice that, differently from the Castellano and 
Pastor-Satorras model [141 ]. the infection spreading occurs among individuals with the same 



number o 
in ref. 



connections (the regular lattice coordination) as in the original CP [16], while 
14| the number of connection varies from one individual to other following a power 
law distribution. Also, the original CP is obtained when a = (no travels), but none limit 
corresponds to the CP on scale-free networks. 

Computer simulations were done for chains of sizes varying from L = 10 2 to 12800 and 
number of nodes ranging from N = 400 to 12800. For each group of 10 samples, a SF 
network was generated from Nq = 10 fully interconnected nodes accordingly the BA prefer- 
ential attachment algorithm jlj. New nodes were sequentially added to the growing network 
through m = 4 links, the minimum degree value. In order to determine the criticality at the 
transition to the absorbing state, the overall density p of infected sites through all nodes, 
the local density ( restricted to those nodes with at least one infected site, the fraction of 
colonized nodes Q (those with infected sites), and the survival probability P s were evaluated. 



III. RESULTS 

As observed for the CP in SF networks jl^|, the simulations in FIG. [2] revealed a non 
null epidemic threshold. The upper inset shows the corresponding survival probabilities, 
exhibiting the same behavior. The initial condition was a single infected site at the center 
of a randomly chosen node. In the upper-critical regime, both quantities reach constant 
asymptotic values, contrasting with the increasing as n ~ t d of the mean number of infected 
sites for the CP on regular lattices [16]. The heterogeneous distributions of infected sites 
in the nodes, with only few of them highly infected, produce the plateaus observed at 
relatively short times (~ 10 3 for chains with L = 6000 sites). In general, models with 
absorbing configurations exhibit power law time dependencies P s ~ t~ s , n ~ t v and R 2 ~ t z 
at the critical point for the survival probability, the mean number of infected sites, and the 

nnn 

spreading of epidemics, respectively [l_6|, [20|, |21]. Accordingly, FIG. [2] provides spreading 
exponents close to zero for the CP model on the present hybrid network. 
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FIG. 2: Evolution in time of the density of infected sites for L = 10 4 , N = 1600, a = 0.5, and 
epidemic rates varying from A = 1.549 to 1.558 in intervals AA = 10 -3 from the bottom to the 
top. Upper inset: corresponding survival probabilities as functions of time. Lower inset: semi-log 
plots for A = 1.549 — 1.551. The averages were done over N s = 10 4 to 5 x 10 5 independent samples 
(the lower A the larger N s ). 



The mean lifetime divergence when A — > A~ can be used to estimate the critical point. 
The exponential decay p ~ exp(— t/r p ) [lower inset of FIG. [2] determines the decay time 
as a function of A. As shown in FIG. [3j an exponential divergence given by t p (A, L) ~ 
exp(const./A) was found. Here, A = |A C — A|. This behavior differs from the power law 
r ~ |A|~ v ll usually observed in transitions to absorbing states in regular lattices [la ]. Indeed, 
a power law also fits the data, but very large exponent values {y\\ ~ 6 to 8) were obtained, 
supporting further the hypothesis of exponential divergence. The exponents v\\ = oo and 
5 = rj = z = are consistent with the well known scaling relations 5 — (3/f\i, z — 2z/j_/z/|| and 
A5 + 2ri = dz [la]. Despite very strong finite size effects, a characteristic feature of epidemics 

n n n 

on SF networks [9|, HU, HJ], the critical rate can be determined by exponential nonlinear 
fittings. In the inset of FIG. El the critical rate is shown function of L-^Q and its 
extrapolation to an infinite size provides a critical rate A c = 1.5432(5), in which the error is 
in parenthesis. This value agrees with estimates from quasi- stationary simulations described 
in the next paragraphs. The critical rate and exponential divergence of t p does not depend 
on network size N for the studied range. The same analysis was done for ( and P s providing 




40 60 80 100 



FIG. 3: Characteristic infection lifetimes in the sub-critical regime against the infection rates for 
several node sizes. N = 1600 and a = 0.5 are fixed. Inset: FSS analysis to extrapolate the critical 
rate. Sampling as in FIG. El 



the same critical rate and exponential divergence for the mean lifetime estimated through 
the analysis of p. 

In order to perform a finite size scaling (FSS) analysis near the criticality, we used qua- 

sistationary (OS) simulations in which only the survival trials are considered in the averages 
I — I I — I 

[14J, H6[. Moreover, a generalized scaling hypothesis that explicitly considers the validity of 



the usual FSS 114, 161 for fixed iV or fixed L was assumed 



p 8 {A, L, N) = L-^iV-W jqAL^iV 1 /^]. (1) 

Similar relations can be proposed for the other quantities. Thus, the thermodynamic limit 
can be estimated through three distinct ways. The first one, in which the size L grows faster 
than the number of nodes N, implying N/L — > 0, seems to be more adequate to describe 
epidemic spreading in the real world. Alternatively, the FSS analysis can be performed with 
N growing faster than L, leading to N/L — > oo, or maintaining a constant ratio N/L. 

For a large but fixed N, Eq. (JTJ becomes p s (A, L) = L" /3/ ^^(AL 1/ ^), in which Q(x) ~ 
x@ for x ^> 1 and Q{x) ~ x~ 1/J - +/3 for x 1 is the scaling function 16j. So, at the critical 
point (A = 0), the density scales as p s ~ L~P/ U± . Similarly, power laws can be associated to 
the other quantities at the criticality, namely, ( s ~ L ' v± , t ~ L u \\/ U± , and Q s ~ L 7 / v± , as 
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FIG. 4: QS densities p s (squares), ( s (circles) and Q s (diamonds), for fixed N = 1600 and a = 0.1 
at the critical rate A c = 2.118. Upper- and a sub-critical data for p are shown. Straight lines are 
power law fits. Lower inset: QS local densities as functions of time. Upper inset: Critical QS 
relaxation time determined using the density Q s . The number of samples varies from N s = 2 x 10 3 
to 10 5 (the larger L the lower N s ). 




shown in FIG. HI Notice that the density of infected nodes fl s grows with the node size. The 



exponents and the critical rates, determined through the null curvature criterion 23], are 
listed in Table HI These values are, within the margins of error, independent on the number 
of nodes for iV = 3200, 6400, and 12800, or travel rates on the studied interval. In turn, 
the critical infection rate decreases continuously from A c = 3.2928 (the CP critical rate) for 
a = to A c = 1 for a — > oo. 

TABLE I: Critical rates and exponents for the CP on hybrid networks with iV = 1600 nodes. The 
uncertainties in the last digits are indicated in parenthesis. 



a A c f3/uj_ 7/W v\\/vx 

0.1 2.1179(5) 0.526 0.913 0.366 0.486 

0.5 1.5435(5) 0.523 0.905 0.365 0.482 

1.0 1.3451(4) 0.523 0.901 0.367 0.488 

Mean value 0.524(3) 0.906(6) 0.366(3) 0.485(3) 



Analogously, for fixed L, Eq. ([!]) leads to similar power laws at the criticality associated 
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to another scaling function ~H(x), except for Q s that now decays as A 7 '/ 1 ^ since 7' < 0. In 
Table HI1 the critical exponents obtained through distinct approaches to the thermodynamic 
limit are compared. These values confirm the critical behavior foreseen by Eq. ([1]) since the 
sum of the exponents obtained for L and for A fixed equals, within the margins of error, 
those with a fixed ratio N/ L. Indeed, 0" = z/|_/3+z/^/3', i/[ = v\_v'^ and so forth for N/ L fixed. 
Also, whatever the approach, the relation # = /3 + 7 is valid, reflecting the independence 
between densities of infected individuals inside the nodes and the fraction of cities in which 
the epidemics persists at criticality. Indeed, as a result, ( s ~ J> s /Q s ~ L~^ +J ^ v± N~^ ,+J ^ u± . 

TABLE II: Critical exponents for the densities of the CP on hybrid networks for a = 0.5 and 
A c = 1.543. Fixed L FSS analysis were done for L = 1600 and 3200 and N > 800. The ratios 
k = 0.5, 1, and 2 were used in the N/L = k approach. 



Approach 




Exponents 


at A = a ) 








0/v ± 


j/u± 


v \\/ v ±- 


fixed N 
fixed L 
fixed N/L 


0.526(3) 
0.912(9) 
1.428(6) 


0.911(6) 
0.073(8) 
0.97(1) 


0.367(3) 
-0.82(1) 
-0.458(8) 


0.484(3) 
0.080(7) 
0.557(6) 



a ) Primes were omitted for sake of brevity. 

Finally, according the scaling hypothesis p]), plots of L /3 ^- L p s (A, L) x AL 1 ^ U± , for fixed 
N, and A"^'/ v J-p s (A, N) x AA 1 /^, for fixed L, using the correct exponents v± or u' ± should 
collapse the data onto universal curves. In FIG. [5] (upper) are shown the collapses of p 
obtained onto Q{x) (fixed N) for l/u± = 0.423(4) and onto %{x) (fixed L) for l/u' ± = 
0.125(5). Thus, using the exponents from Table HT1 we found (3 = 1.23(3), 7 = 0.86(3) and 

= 1.14(3) for fixed N, and 0' = 7.3(3), 7' = 0.6(3) and vt = 0.6(3) for fixed L. The large 
value found for representing a rapid vanishing of p at Apis consistent with the exponential 
approach obtained for the SIS model in BA networks 9), Also, the finite value obtained 
for Pn in the QS analysis for fixed N deserves a comment. It contrasts with the exponential 
divergence of the characteristic epidemic mean lifetime (FIG. [3]). This paradoxical result 
could be interpreted in terms of two time scales for the infection dynamics. The first one 
refers to an exponentially long persistence time of the epidemics spreading from a spatially 
concentrated focus in the subcritical phase. We conjecture that the heterogeneity of the 
connection among cities, thereby of travels, promotes this long persistence, in contrast to 



10 



10 



10 



10 = 



10 



• N=800 
■ N=1600 

♦ N=3200 
i N=6400 
o L=100 

□ L=200 

o L=400 

a L=80O 

< L=1600 



10 



10 



10 



10 



Q. 10 z 



10 



• N=800 
■ N=1600 

♦ N=3200 
» N=6400 
□ L=100 

L=200 
A L=400 
< L=800 

1 L=160O 



1A» 



3'* 



V 



10 



10 



10 



FIG. 5: Upper: collapses of p for either N = 1600 or L = 800 fixed. For fixed N, l/u± = 0.423 
and /3/v± taken from Table ITU were used. For L fixed, 1/V^_ = 0.125 was used. Bottom: collapses 
for both L and N varying, a = 0.1 in this Fig.. Here, p* = p s L^/ v± or p s N^'^ and A* = AL 1 /^ 
or AJV 1 /"!. In turn, p** = p s L^l v± N^/ v 'x- and A** = AL^-L^K 



the CP with sources 



16j if the SF network is replaced by a regular one. The second time 



scale is associated to the outbreaks, which diverges algebraically at the critical point as in 
local epidemics. The exponents listed in Tables 1 and 2 are quite distinct from those related 
to DP universality class (fixed N) and, also, do not agree with numerical estimates for 
CP on SF networks l^j (fixed L), establishing a new universality class. 
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IV. SUMMARY 



We presented numerical simulations of the CP on hybrid networks in which the nodes 
are themselves regular lattices and individuals travel between nodes. The main results are a 
finite epidemic threshold and an exponentially long persistence of the epidemics below the 
critical rate, concomitantly with a power law divergence of the outbreak's duration. Gener- 
alized scaling functions involving both regular and SF components were proposed and the 
associated critical exponents determined, demonstrating that the CP on this hybrid net- 
work and nonvanishing travel rates establishes a new universality class. Also, at criticality, 
the densities of infected individuals inside the nodes and the fraction of cities in which the 
epidemics persists are independent. 
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